Dimensional Crossover of Dilute Neon inside Infinitely Long 
Single- Walled Carbon Nanotubes Viewed from Specific Heats 

Zhan-chun Tu 1,2, Q and Zhong-can Ou-Yang 1 ' 3 

1 Institute of Theoretical Physics, The Chinese Academy of Sciences, 

P.O.Box 2735 Beijing 100080, China 

2 Graduate School, The Chinese Academy of Sciences, Beijing, China 

3 Center for Advanced Study, Tsinghua University, Beijing 100084, China 

Abstract 

A simple formula for coordinates of carbon atoms in a unit cell of a single-walled nanotube 
(SWNT) is presented and the potential of neon (Ne) inside an infinitely long SWNT is analytically 
derived under the assumption of pair- wise Lennard- Jones potential between Ne and carbon atoms. 
Specific heats of dilute Ne inside infinitely long (5, 5), (10, 10), (15, 15) and (20, 20) SWNT's are 
calculated at different temperatures. It is found that Ne inside four kinds of nanotubes exhibits 
3-dimensional (3D) gas behavior at high temperature but different behaviors at low temperature: 
Ne inside (5, 5) nanotube behaves as ID gas but inside (10, 10), (15, 15), and (20, 20) nanotubes 
behaves as 2D gas. Furthermore, at ultra low temperature, Ne inside (5, 5) nanotube still displays 
ID behavior but inside (10, 10), (15, 15), and (20, 20) nanotubes behaves as lattice gas. 

PACS numbers: 61.46. +w, 82.60.Fa 



Since the discovery of carbon nanotubes,— the peculiar electronic and mechanical prop- 
erties of these structures have attracted much attention^M Experiments have also revealed 
that carbon nanotubes can be used to store hydrogen^ and other gases.— Many physicists 
expect gases in nanotubes or nanotube bundles to display 1-dimensional (ID) behavior as 
a consequence of the remarkable aspect ratio of the length of tubes to their radius. Cole 
and his co-workers^ and other researchers^ have theoretically studied properties of gases in 
nanotubes or nanotube bundles. One of the most fantastic properties they found is specific 
heat of dilute gas inside single-walled carbon nanotubes (SWNT's) as a functions of temper- 
ature: With temperature increasing it shows the thermal behavior changing from ID to 2D, 
to 3D. However, there is still a question: Do ID and 2D behaviors always exist for dilute 
gas inside SWNT's if the temperature is low enough? 

A series of skillful experiments to measure specific heat of gas in carbon nanotube bundles 
were done by Lasjaunias et a?.— In their experiments, He atoms are adsorbed within the 
interstitials or external grooves and surfaces of the bundles, and exhibit ID behavior at low 
temperature. Their experiments suggest that the specific heat can reflect the dimensional 
information well. 

Here we firstly present a simple formula for coordinates of carbon atoms in a unit cell of a 
SWNT and derive an analytical expression of the potential of neon (Ne) inside an infinitely 
long SWNT under the assumption of pair-wise Lennard- Jones potential between Ne and 
carbon atoms. And then, we calculate specific heats of dilute Ne inside (5, 5), (10, 10), (15, 
15) and (20, 20) SWNT's and find that Ne inside four kinds of nanotubes exhibits 3D gas 
behavior at high temperature but different behaviors at low temperature: Ne inside (5, 5) 
nanotube behaves as ID gas but inside (10, 10), (15, 15), and (20, 20) nanotubes behaves 
as 2D gas. Furthermore, at ultra low temperature, Ne inside (5, 5) nanotube still displays 
ID behavior but inside (10, 10), (15, 15), and (20, 20) nanotubes behaves as lattice gas. Ne 
inside (5, 5) nanotube does not display 2D behavior and inside (10, 10), (15, 15), and (20, 
20) nanotubes does not display ID behavior. Thus Ne inside nanotubes does not always 
display ID behavior in spite of the remarkable aspect ratio of the length of tubes to their 
radius. Besides, Ne inside a definite nanotube may display neither ID nor 2D behaviors at 
different temperatures. 

A SWNT without two caps can be constructed by wrapping up a single sheet of graphite 
such that two equivalent sites of hexagonal lattice coincided To describe the SWNT, some 
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characteristic vectors require introducing. As shown in Fig^ the chiral vector C^, which 
defines the relative location of two sites, is specified by a pair of integers (ni,ri2) which is 
called the index of the SWNT and relates Ch to two unit vectors &i and a 2 of graphite 
(Ch = niB-i + n 2 a 2 ). The translational vector T corresponds to the first lattice point of 
2D graphitic sheet through which the line normal to the chiral vector Ch passes. The unit 
cell of the SWNT is the rectangle defined by vectors Ch and T, while vectors ai and a 2 
define the area of the unit cell of 2D graphite. The number N of hexagons per unit cell of 
SWNT is obtained as a function of n\ and n 2 as N = 2{n\ + n\ + r^n^/fiR, where is the 
greatest common divisor of (2n 2 + rii) and {2n\ + n 2 ). There are 2N carbon atoms in each 
unit cell of SWNT because every hexagon contains two atoms. To denote the 2N atoms, 
we use a symmetry vector R to generate coordinates of carbon atoms in the nanotube and 
define it as the site vector having the smallest component in the direction of Ch- From a 
geometric standpoint, vector R consists of a rotation around the nanotube axis by an angle 
= 2tt/N combined with a translation r in the direction of T; therefore, R can be denoted 
by R = (^(t). Using the symmetry vector R, we can divide the 2N carbon atoms in the 
unit cell of SWNT into two groups: one includes N atoms whose site vectors satisfy 

Ai = IK — [IH ■ T/T 2 ]T (I = 0,1,2, •• •, N- 1), (1) 

another includes the remainder iV atoms whose site vectors satisfy 

B ; = /R + B - [(ZR + B ) ■ T/T 2 ]T 

- [(/R + B ) • Ch/C 2 h }C h (I = 0, 1, 2, • • • , N - 1), (2) 

where Bo represents one of the nearest neighbor atoms to Ao- In and only in above two 
equations, [• • •] denotes the Gaussian function, e.g., [5.3] = 5. 

To obtain the potential of Ne inside the nanotube, we firstly consider another simple 
system shown in Fig. El Many carbon atoms distributed regularly in a line form an infinite 
atom chain and a Ne atom Q is out of the chain. The interval between neighbor atoms in the 
chain is T, and the site of atom Q relative to atom can be represented by two numbers c% 
and c 2 . We take the Lennard- Jones potential U(Rj) = 4e[(cr / Rj) 12 — (cr/Rj) 6 } between atom 
Q and atom j in the chain, where Rj is the distance between Q and atom j, and e = y/e c e ne , 
a = (a c + a ne )/2 with e ne = 35.6 K, a ne = 2.75 A, e c = 28 K and a c = 3.4 A.^ii We express 



3 



the potential between atom Q and the chain as 

U QO = 4e[a 12 f/ 6 ( Cl , c 2 ) - a 6 U 3 (c u c 2 )], (3) 

oo 

where £4(ci,c 2 ) = Yl [( Cl +nT) 2 +c i ] k = lj 2, • • •) which can be calculated through the 

n=— oo 2 

following recursion^ 2 - 

TT (r r\ — 7rsinh(27rc 2 /T) 

vi\yi,L,2) C2 T[cosh(27rC2/r)-cos(27rci/T)] ' 

U k+1 (ci, c 2 ) = -l/(2kc 2 )dU k /dc 2 . 
The (20, 20) tube, for example, with infinite length can be regarded as 2N = 80 chains. 
Thus the potential of any point Q inside the tube can be calculated as 

2N 

U(r,e,z)=J2u QC , (5) 
i=i 

where (r, 9, z) is coordinates of Q in the cylindrical coordinate system whose z-axis is the 
tube axis, r the distance between Q and z-axis, and 9 the angle rotating around z-axis from 
an axis which is vertical to z-axis and passes through atom A on the tube to the plane 
that contains Q and z-axis. As an approximation, we neglect the potential varying with z 
and 9 because we find that it is much smaller than the potential varying with r through 
calculations, and fit the potential with U(r) = Ae[(-^) 10 — (-^) 4 ), where p = 13.56 A is 
the radius of the tube, e = 390 K, and a = 2.63 A (see also FigJBJ). Moreover, we simplify 
it as 

f 4e[(_^)io _ (_2_)4i r<p -v 
U(r) = I [Kp - r) Kp - r) 1 H ' (6) 

I oo, r > p — a. 

Because we consider the dilute Ne, we can neglect the interaction between Ne atoms and 
write the single particle Schrodinger equation^ as Hip = Eijj, where H = — f^V 2 + U(r) 
and ip = ( f )e l ( me + KZ ) . From them we arrive at 

E = ^ + E m (KGR,m = 0,±l,±2,...), 

H{r)<t> = E m <f) } (7) 



Setting r = (p — er)£, e = and r] = a/p, Eqs.(jnj) and (JJJ) are transformed into 



u(0={ ^-V-VS' ^-V-vK ' J) ' ' ( 8 ) 
00, £>1. 



and 

If we let \<f) = J2 n a n\Xn) , we will obtain the secular equation 

det(H jn - E rn S jn ) = 0, (10) 

where H jn = J Xj^iOXn^, and S jn = XjXn£d£- If we let Xn = J\m\(v n 0, where J\m\(0 
is the m-th order Bessel function of the first class and v n is the n-th zero point of Bessel 
function,— we can calculated E mn , (m = 0, ±1, ±2, • • • ; n = 1, 2, 3, . . .) from Eq.([TU|). 

If there are M Ne atoms inside the tube, we have the free energy T = —AfTlnZ, where 

ft 2 K 2 

T is the temperature, Z = e~ Emn ^ T e~ 2 ^ r dn and the Boltzmann factor is set to 
l.J£ We can easily obtain the specific heat per atom is 



MdT 2 2 T 2 

where (£> = E ^7- e EW T and (E 2 ) = . 

In Fig0J the symbols "V" reflect c„ varying with temperature T, which implies Ne inside 
(20,20) tube behaves as 3D gas at high temperature (specific heat approaches to 3/2) and 
2D gas at low temperature (specific heat is 1). Therefore, we can naturally assume that all 
atoms in the valley of potential U(r) at low temperature, i.e. lie on the shell S* with radius 
f? = p[l-(5/2) 1 /6 r? ]. 

Now we consider the thermal property of Ne inside the nanotube at low temperature more 
carefully. From Eqs.(j2J)-© we can easily calculate the potential U s (v, z) on S*, where v = q6. 
The Hamiltonian of single particle can be expressed as H' = — f^(J^ + J-js) + U 8 (v, z). 

In fact, U s has periodic structure. If we denote Qi = (2ttq/N,t), ol 2 = (0,T) and 
ji = 1\oli + l 2 et 2 , we know U s (r + 7^) = U s (r), where li, l 2 G Z and r = (v, z). On the one 
hand, we have the Bloch's theorem*^ 

r) r), 
$(k, r + ji) = e iK ' 7i $(n, r), 

which suggests that $(k, r) = a(K+Gj)e^ K+GjS) ' r , where Gj = ji{3i+j 2 {3 2 with ji, j 2 G Z, 
f3\ = (N/g,0) and f3 2 = (—tN/(Tq),2tt/T). From Eq. (fT2"j) . we obtain the secular equation 



det{H l3 - EJu) = 0, (13) 
5 



where Uij = ^{k + G 3 -) 2 5zj + U tj , U tj = ^ f Qo e^ Gl - G ^ r U s {r)dr and fi = |«i X a 2 \. 

On the other hand, Periodic boundary condition along the circumference of the shell 
S* suggests that we just need to consider the first Brillouin zone which consists of n = 
(m„/ g, k z ) where m v G Z, < m v < N and k z G R, < k z < 2%/T. 

From Eq. (jl3|) . we can calculate the energy E mvjKz , and then the free energy T = 
—NThiZ?L where Z = E m e~ Emv ' K - z ^ r dn z . Moreover, the specific heat per atom 

is±£ 

v ~ Afdr* ~ r 2 ' [ ' 

TV— 1 N—l 

where (E') = E E mv ^e~ E ^^ ' r dK z and (E> 2 ) = E fi*' T < >Kz ^— /T ^- 

m v =0 m v =0 

In FigEl the symbols "V" reflect c v varying with the temperature T, which implies Ne 
atoms inside (20,20) tube behave as the lattice gasA& at ultra low temperature (specific heat 
is 0) and 2D gas at low temperature (specific heat approaches 1). There is no ID gas inside 
(20,20) tube, which is quite different from our usual notion. 

Adopting the similar method, we calculate c v of Ne inside (5,5), (10,10) and (15,15) 
nanotubes and show them in Fig0] and Fig|SJ When we obtain FigJU we have adopted the 
approximate potentials as U(r) = 4e[(^) Pl — {-^) P2 }. The parameters e, a, p, pi, and 
p 2 are shown in Table U] These two figures imply that Ne inside four kinds of nanotubes 
exhibits 3D gas behavior at high temperature but different behaviors at low temperature: 
Ne inside (5, 5) nanotube behaves as ID gas (specific heat is 1/2) but inside (10, 10), (15, 
15), and (20, 20) nanotubes behaves as 2D gas. Furthermore, at ultra low temperature, Ne 
inside (5, 5) nanotube still displays ID behavior but inside (10, 10), (15, 15), and (20, 20) 
nanotubes behaves as lattice gas. Ne inside (5, 5) nanotube does not display 2D behavior 
and inside (10, 10), (15, 15), and (20, 20) nanotubes does not display ID behavior. Thus 
Ne inside nanotubes does not always display ID behavior in spite of the remarkable aspect 
ratio of the length of tubes to their radius. Besides, Ne inside a definite nanotube may not 
display both ID and 2D behaviors at different temperatures. 

Our calculations support the view of dimensional crossover of gas inside SWNT's put 
forward by Cole et.al^ but there is a significant difference: Cole et.al. obtain the dimension- 
ality varying from one to two to three with the temperature increasing, but we obtain the 
dimensionality varying either from one directly to three or from two to three. The reason is 
that the potential of gas inside SWNT's selected by Cole et.al. departs father from the real 
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case than that done by us. 

If the potential on shell S* is uniform, one can intuitively expect the dimensionality to 
vary from one to two to three as the axial, azimuthal and radial degrees of freedom become 
excited as temperature increases. But the potential on this shell has more exquisite structure 
(the periodic structure) for large diameter nanotubes, such as (10,10), (15,15) and (20,20) 
tubes, which results in the approximate same axial and azimuthal excited energy. Therefore 
Ne in large diameter nanotubes exhibits lattice gas (0D), 2D and 3D behavior as temperature 
increases. On the contrary, the shell converges to the tube axis for small nanotube, e.g. (5,5) 
tube, which results in the dimensionality varying from one direct to three as temperature 
increases. 

The authors acknowledge the useful discussions with Prof. H. W. Peng, Dr. M. Li, F. 
Liu, Y. Zhang, L. R. Dai, and J. J. Zhou. Furthermore, we thank Dr. H. J. Zhou because 
he help us read the English. 
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FIG. 1: The unrolled honeycomb lattice of a SWNT. By rolling up the sheet along the chiral 
vector C/j, that is, such that the point Aq coincides with the point corresponding to vector C^, a 
nanotube is formed. The vectors ai and SL2 are the real space unit vectors of the hexagonal lattice. 
The translational vector T is perpendicular to Ch and runs in the direction of the tube axis. The 
vector R is the symmetry vector. Aq, Bq and Ai,Bi(l = 1, 2, • • • , N) are used to denote the sites 
of carbon atoms. 
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FIG. 2: An infinite atom chain and an atom Q out of the chain. Many atoms distribute regularly 
in a line form the infinite atom chain. The interval between neighbor atoms in the chain is T, and 
the site of atom Q relative to atom can be represented by numbers c\ and C2- 



TABLE I: The parameters in the approximate potentials for different nanotubes. 

(n,m) e (K) a (A) p (A) Pl p 2 

(5,5) 520 2.63 3.39 10 1 

(10,10) 360 2.63 6.78 10 3 

(15,15) 418 2.64 10.17 10 4 
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FIG. 3: The potentials inside the (20,20) nanotube calculated from Eq.@ (squires) and the fit 
curve (solid curve). 
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FIG. 4: The specific heats per atom c v of Ne inside (5, 5), (10, 10), (15, 15) and (20, 20) nanotubes 
under the approximate potentials varying with the temperature. 
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FIG. 5: The specific heat per atom c v of Ne inside (5, 5), (10, 10), (15, 15) and (20, 20) nanotubes 
at low temperatures. 
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